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ABSTRACT 



The technique of model atmosphere calculation for magnetic Ap and Bp stars with polarized radiative transfer and magnetic line blanketing is 
^SJ ■ presented. A grid of model atmospheres of A and B stars are computed. These calculations are based on direct treatment of the opacities due to 
^ ■ the bound-bound transitions that ensures an accurate and detailed description of the line absorption and anomalous Zeeman splitting. The set 
] of model atmospheres was calculated for the field strengths between 1 and 40 kG. The high-resolution energy distribution, photometric colors 
\l and the hydrogen Balmer line profiles are computed for magnetic stars with different metallicities and are compared to those of non-magnetic 
reference models and to the previous paper of this series. The results of modelling confirmed the main outcomes of the previous study: energy 
] , redistribution from UV to the visual region and flux depression at 5200 A. However, we found that effects of enhanced line blanketing when 
transfer for polarized radiation takes place are smaller in comparison to those obtained in our first paper where polarized radiative transfer was 
neglected. Also we found that the peculiar photometric parameter Aa is not able to clearly distinguish stellar atmospheres with abundances other 
than solar, and is less sensitive than A{Vi - G) or Z to a magnetic field for low effective temperature (retf = 8000 K). Moreover we found that 
the back determination of the fundamental stellar atmosphere parameters using synthetic Stromgren photometry does not result in significant 
errors. 
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J> 1. Introduction Imaging technique (including all Stokes parameters). An enor- 
mous part of this work is currently being performed using clas- 

The magnetic line blanketing in stellar atmospheres due to the . j i , . 

^ SIC model atmospheres. 

_j , Zeeman eiiect is supposed to be responsible for observable . , , • ^ , , ^ 

, ■ . • , • ,. /^T,N In general, the magnetic held afreets the energy transport, 
- - ' characteristic features of magnetic chemically peculiar (CP) 

stars, for instance flux depressions (Kodaira 1969) in the vi- 
sual sp ectrum, and flux redistribution from UV region to the ^ , , ^ ... . 

1 — — » J- 1- ^1 rect, due to complex Zeeman splitting, assuming a simplihed 

visual (Leckrone 19731. A review of hterature of early at- ^ , ^^^.^ ^ ,^ .„ ^;„„i„,°„ tu„ i; f 

tempts to consider effects of magnetic field in atmospheres 
modelling, as well as a simplified but state-of-the-art technique 
of stellar model atmospheres calculation taking into account 



hydrostatic equilibrium, diffusion processes and the formation 
of spectral lines. In Paper I we considered in detail the last ef- 



model of a magnetic field vector perpendicular to the line of 
sight and estimated its influence on model atmospheres with- 
out solving the polarized radiation equation. 



u » J • • 1 This paper continues the investigation and aim started in 

the magnetic absorption, were presented i n our previous works ^ ^ ° 

JKhan et alJ2004HKochukhov et alhool hereafter Paper I). *^ P^P^' °f ^^'^^^^^^ ^ realistic calculation 

™,, .r- jji. 1. £ r^n of the anomalous Zeeman effect in the classical 1-D models 

The development of improved model atmospheres for CP „ 

, . , , , 1 r J . J- of stellar atmospheres and to investigate the resulting effects 

stars IS an extremely important task for understanding the na- ^ o o 

c c ^ c i-^T, » i/rji» u on the model structure, energy distribution and other general 

ture of features of CP stars. Model atmospheres are used for ° 

■ , n ^ „ ^ . • ,1 ■ V ^- ^ observables. In this work we summarize the transfer equation 

a variety of stellar astrophysics problems in application to CP ^ 

^ , .j^ -.- 1.J 1 for polarized radiation, describe the general technique of our 

stars: fundamental parameters determination, abundance analy- ^ ' o -i 

^ „ ^- £ ij . 1. J ^ 1 J 1- computational procedure, and calculate the same grid of model 

SIS, stellar magnetic field geometry research, detailed fine pro- ±- ±- ' o 

, J /• 1 J- r 11 . . . £ ^i. rr a: .\ atmosphercs of A and B stars as in Paper I exploring a Tgfi 
file study (including the full treatment of the Zeeman effect), ^ i- ±- o cu 
, ^ „ £ ^ ^- 1. .1, 1 range relevant for CP stars. The paper is organized as follows, 
and stellar surface properties reconstruction by the Doppler ^ m , . ,-7, 
In Sect.|2|we describe our technique of model atmosphere cal- 

Send ojfprint requests to: S. A. Khan, e-mail: skhanSastro . uwo . ca culation and numerical implementation of the transfer equation 
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for polarized radiation in the line opacity calculation. Section|3] 
presents results of numerical application and, Sect.|4] summa- 
rizes our work. 



2. Calculation of magnetic model atmospheres 

The model atmospheres calculated in this paper were con- 
structed using a m odified version of the L LModels code, origi- 
nally developed bv LShulvak et al.l ll20Q4l) . The LLModels code 
is an LTE (local thermodynamic al equilibrium) 1-D stellar 
model atmosphere code for early and intermediate type stars 
which is intended for as accurate a treatment as possible of 
the line opacity using a direct method for the line blanketing 
calculation. This approach allows us to take into account indi- 
vidual chemical compositions (not scaled to the solar) of stel- 
lar atmospheres, possibly inhomogeneous vertical distribution 
of abundances (as a result of diffusion processes which we do 
not consider here in a self-consistent way), and individual line 
broadening mechanisms such as anomalous Zeeman splitting 
of spectral lines that are required for proper modelling of at- 
mospheres of CP stars. 

Using other techniques (Opacity Distribution Function, 
Opacity Sampling) for modelling atmospheres of CP stars is 
very complicated or even impossible due to their basic limi- 
tations and statistical nature. These methods were developed 
several decades ago in order to reduce the computational time 
of the treatment of millions of spectral lines. 

In contrast, the direct method implemented in the 
LLModels code is free of any approximations so that it fully 
describes the dependence of line absorption coefficient upon 
frequencies and depths in a model atmosphere, it does not re- 
quire pre-calculated opacity tables, it has no limits on number 
of frequency points or number of spectral lines, and it has very 
reasonable execution time on an ordinary PC. 

The c ode is based on modified ATLAS9 subroutines 
llKuruczl l993a) and on the spectrum synthesis code described 
bylTsvmbal ( 1996). In the new version of the LLModels code, 
the continuum opacity sources and partition fu nctions of iron- 
peak elements from ATLAS 12 llKuruczll993hl) ar e used. Parts 
of the magnetic spectrum synthesis code SynthM (l Kharj2004 
were incorporated into LLModels to deal with transfer of po- 
larized radiation. 



2.1. Transfer equation for polarized radiation 

In this subsection we summarize the basics of radiative trans- 
fer in the presence of the magnetic field, formulated in suit- 
able way for our calculations, and point out some of its im- 
portant details. The comprehensive description and math for 
polarized radiation can be found in several books (Rees 19 87; 
iRees & Mumhvll987tlst^flnll994l: ldd" Toro Iniestal2nn3h . 

The polarized beam can be fully described in term of Stokes 
formalism 



/ = (/, e, u, vy 



(1) 



where / is called Stokes vector. The parameter / represents the 
total intensity of the beam, Q and U the difference between 



the intensities of linearly polarized components along position 
angles (Q,n/2) and {n/4,3n/4), and V the intensity difference 
between right-handed and left-handed circular polarized com- 
ponents. 

The Stokes formulation of the transfer equation on a scale 
of column mass m (which is usually used in model atmospheres 
calculation) is 



am 



(2) 



where yu = cos(6') is the cosine of the angle 6 between the light 
propagation direction and the perpendicular to the surface, K is 
the total propagation matrix for all contributing lines. 



(3) 



If we assume unpolarized continuum radiation and consider co- 
herent isotropic scattering, neglecting polarization effects, then 
the total emission vector j is 



j = /fc 5c lo + 5l1o ^ fo""^!™ + 0-cJvlo- 



(4) 



lines 



Here 1 is a 4 x 4 identity matrix, Iq - (1,0, 0, 0) ' , and are 
the continuum mass absorption coefficient and source function, 
ATo is the line center mass absorption coefficient for zero damp- 
ing, zero magnetic field and corrected for stimulated emission, 
cr^ is continuous scattering coefficient, and Jy = ^lydco/An is 
mean intensity. 

If we assume LTE then, - Sl = By, the Planck function 
at the local electron temperature. In that case introducing the 
source vector 



5 = (1 - q;)Bv1o + Q' Jylo, 

where a - cTcIki, we can rewrite Eq. as 

am 



(5) 



(6) 



To comput e the source ve ctor S we perform A-iterations as de- 
scr ibed bv iKuruczllll97(]h . 

Finally, the line propagation matrix O is 



' 0/ 4>Q 


4>u 


<Av ^ 




ipv 














4>i ' 



(7) 



where 






<A/ = 


5(0/. 


sin^ 7 + :j(0r + <Ai)(l + cos^ 


<PQ = 




- 5(0,. + (pb)) sin^ r cos 2x, 


4>V = 


Wp 


- j((/>r + (f>b)) sin^r sin2;t'. 


H - 




- (I>b) cos y, 


4lQ = 


5(<Ap 


- 5(1/',. + if/b)) sin^ r cos 2x, 


lAt/ = 




-\{^r + 4fb)) sin^ r sin 2x, 


lAv = 


5('A. 


- i/'fc) cos y. 



(8) 
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Now, let us return to the absorption 4>p,b,r and the anomalous 
dispersion 4'p.h.r profile definitions. In the presence of a mag- 
netic field an atomic level k defined by the quantum numbers 
Jk, Lk, S It splits into the 27^ + 1 states with the magnetic quan- 
tum numbers - -Jk, . . . , +Jk- The absolute value of the 
splitting is defined by the field modulus |B| and by the Lande 
factor ^t, which in the case of LS coupling can be calculated as 



3 Sk{Sk + l)-U{U + l) 



2Jk(Jk + 1) 



(12) 



Fig. 1. Reference frame for representing the Stokes vector / 
and magnetic field vector B. The z-axis shows the direction of 
the beam propagation (or line of sight). The magnetic vector 
has an inclination y to the propagation direction and azimuth 
measured counter-clockwise from the x-axis: < y < tt, < 
X ^ 27r. 

The angles y andx are defined in the cartesian reference frame 
xyz as shown in Fig.^ The <pp,h.r and ^p^b,r are absorption and 
anomalous dispersion profiles that will be described below. 

For a better understanding of the physical meaning of the 
O matrix, it can be presented as sum of three matrices: 



According to the selection rules the following transitions are 
allowed between the split upper u and lower I levels 



AM = M„ - M/ = 




(13) 



The wavelength shift of the component ij (where /' - p, b, r, 
and ip^b,r - 1, . ■ . ,Np^b,r\ Np^b,r is the number of components in 
each group) relative to the laboratory line centre Ao is defined 
by 







eAl\B\ 
Anmc^ 



(giMi - guMu)i^ 



(14) 



/1 + 



4>Q cf>u 4>v ' 

0e 

0t/ 





(0 Q ^ 

lAv -^u 

-lAi/ lAe 

(Ac -lAe 



(9) 



The relative strengths S ij of the n and cr components are 
given bv lSobelmanI ( 1 1 9771) (see Table 2 in Paper I) and are pro- 
portional to 



/„ 1 Ji 
-M„ M„ - M, Ml 



(15) 



Np Ni, N, 



(=1 



1=1 



The first, diagonal, matrix corresponds to absorption by the 
medium that results in equal changes of all Stokes parame- 
ters. The second, symmetric, matrix corresponds to dichroism, 
where some polarized components are extinguished more than 
others. The third, antisymmetric matrix corresponds to disper- 
sion phenomena (or birefringence), where phase shifts change 
states of linear polarization among themselves (Faraday rota- 
tion) and states of linear polarization with states of circular po- 
larization (Faraday pulsation). The sum of two first matrices 
represents absorption effects, the last matrix represents anoma- 
lous dispersion (or magneto-optical effects). 

The choice of x and y axes in the reference frame xyz 
(Fig.n* is free. In fact the intensity of the beam and degrees 
of linear + U^II and of circular VII polarization are in- 
dependent of the choice of x and y axes. Specifically, the ro- n 
tation of the reference frame by an angle a measured in same ^ = 2 5 , F{a v — vi) 
direction as x modifies the Stokes vector / as i i i 



a 3j - symbol. The normalization of the relative strengths is 
performed to unity for each group of the Zeeman components 



(16) 



The absorption and anomalous dispersion profiles in absence 
of macroscopic velocity field are 



(17) 



and 



i,=i 



(18) 



/' = GL 



where 



G(a) = 



1 0^ 

cos2q' sin2Q' 

- sin 2a cos 2a 

1 



(10) 



(11) 



where H{ a, u) and F(a , u) are the Voigt and Faraday-Voigt 
functions (IStenfloll994 p. 65) 



a r e-y 

Hia,v) = - 

^ J (v-yY + < 



■dy 



(19) 



Hence neither the total beam intensity I nor the V parameter 
depend on the xy choice, whereas Q' and U' are defined by 
simple linear transformations. 



and 



F(a, v) ■■ 



+ 00 

2nj (u- 



y)e 



yY 



■Ay. 



(20) 
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The parameters a, v and u, are expressed in units of the 
Doppler width A/Id 



a = 



(21) 



where T is the total line damping parameter, 
v^iA- Ao)/AAu, 
where A is the current wavelength in the line, and 
y,- = AAiJAAo. (23) 

2.2. Solution of the transfer equation 

Usually, one of two powerful methods to solve polarized trans- 
fer equation are used: the integral Feautrier method and the dif- 
ferential DELO (Diagonal Element Lambda Operator) method 
( Auer et al. 1977; Rees et al. 1989; Piskunov & Kochukhov 
l2002h . which have proved their capabilities for accurate 
solution of the transfer equation. Howev er, the Feautrier 
metho d is only half as efficie nt as DELO ( lRees_eLaD llSSSt 
IPiskunov & Kochukhovlf2002h . Considering this fact, and the 
requirement that for stellar atmosphere modelling in compar- 
ison to spectrum synthesis, one needs to know not only the 
emergent Stokes vector at the surface but the radiation field 
(flux and mean intensity) through out the atmosphere, that re- 
quires the solution of an additional differential equation (and 
extra computational time), the method of choice is DELO. 

The DELO method is a one-way method that allows one 
to calculate inward and outward radiation separately. This 
method, as a differential method, does not require special prop- 
erties of the propagation matrix: the presence of symmetry 
in dire ctions and frequencies. It was developed by Rees et al] 
( Il989l) who used the fact that the diagonal elements of matrix 
K are equal to 



Kl 



lines 



line j^line 

fi ■ 



Defining a vertical optical depth as 
dr = K; dm. 



(24) 



(25) 



a modified propagation matrix K' with zeros on its diagonal 
K' = K/k, - 1 (26) 
and a modified vector 



(27) 



which has the units of a source function, we can rewrite Eq. Q 
in the form 



where 

^^S'-K'I 



(28) 



(29) 



is the new source term (or effective source function). 

Let us introduce a grid of depth points from the surface to 
the deepest layer t, (i - 1, ...,A^) along the vertical direction. 



The formal integral solution of Eq. ( I28> can be written as rela- 
tion between the Stokes vectors /, and 



Ti+l 

/,■ = e,/,+i + J e-(^-^')^(T)dT, 



^■^^^ using a simple trapezoidal formula 



6i = 5(a:/+' + K\){mi+i - mi)lii. 



(30) 



(31) 



In order to solve the integral ( I30> analytically two ap- 
pr oaches were devel oped. The first one, originally proposed 
by iRees et all lll989l) . is based on linear interpolation of the 
source term ^ between grid points, while the second one relies 
on linear i nterpolation only for K'l a nd parabolic approxima- 
tion for 5' (i ocas-Navarro et al .I2OOOI) . The solution of Eq. ( I30> 
for these two cases can be written in common form 



^ili - Yili+i + Zi, 

where 



(32) 



(33) 



and 

ff; = 1 - e,-, 

lii^[\-{\+6dei\l6i. 

For a linear approximation of the source term S' 
Zi = (a,-yS,)5;+AS;^l, 
and for a quadratic approximation 

Zi = a,S;_i + biS'i + c,S-+i, (36) 
where (see lPiskunov & Kochukhovl2002h 



(34) 



(35) 



z - 6i+iy 
(5;+i +6i)y-z 



Ci - x + 



z - (6i + 26i+i)y 



(37) 



X - I - £i, 

y = 5i- X, 
z^6]-ly. 

We used the quadratic DELO meth od because it is significantly 
more accurate that linear DE LO jSocas-Navarro et aLlbOQOl: 
IPiskunov & Kochukhov'2002V 

The recursive Eq. J32> is used to calculate the Stokes vec- 
tors of radiation directed outwards or inwards I through the 
atmosphere with different boundary conditions. To calculate 
we start at the lower boundary assuming unpolarized thermal- 
ized radi ation 1\, - Bv ^l o, or using the improved asymptotic 
formula jAuer et alJl977h 



+ -1, (dSv 



(38) 
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and iterating outwards to the surface. The calculation of the / 
starts with initial estimate 



r(0) = 0, 



(39) 



assuming that there is no inward radiation at the surface, and 
iterates inwards to the atmosphere. 

2.3. Model convergence criteria 

To control convergence of the model calculations we use two 
criteria: the c ondition of rad iative equilibrium a nd consta ncy of 
the total flux dTsvmbal & Sh ulvak 2003; Shul vaket al.ir2 004). 
Both criteria are checked on each iteration and for each atmo- 
spheric layer, and the iterative procedure continues as long as 
errors are more then 1 % and and temperature correction exceed 
IK. 

The condition of the radiative equilibrium (or total energy 
balance) expressing the absence of sources and sinks of energy 
may be easily derived from Eq. (|6jl by integration over solid 
angle and frequency. 



oo 

dm J 



CO 



Ky(Iy -Sy) dCO dV. 



ly COS 6d(L)dv ■■ 



The left side of the first equation of this system is 



oo 

dm J ^ 



ly COS fldw dv - - — F(m) - 0, 
dm 



where F{m) is the total flux at depth m in an atmosphere. Thus 
the right side of the first equation of the system J40> leads to 
formulation of radiative equilibrium as 

CO CO 

^ §(kiI + kqQ + KuU + KvV) du) dv - ^ j> ki Sydcjdv. (42) 



Here we note that the integrand expression of Eq. (142 > does 
not depend on reference frame rotation, which followes from 
Eqs. (0 and O- 

The flux conservation condition is formulated in the usual 
manner 



oo 

/ 



fvdv-o-r^jj. = 0, 



(43) 



where Fy can be calculated as an intensity integral over solid 
angle dcj - smOdOdip - -djj dip as follows 



2n +1 2jt +1 

J J jdl^dfidp - J" J jdl dfidif. 



(44) 



where and I are the first elements of corresponding Stokes 
vectors /■*" and /" which are calculated in Eq. i32\ with bound- 
ary conditions ( I38> and (13 9> . and angles 6 and tp are represented 
in Fig. 121 Here we neglect convective flux because the strong 
global magnetic field may be supposed to prevent turbulent mo- 
tions in the atmospheres of Ap and Bp stars. 




Fig. 2. Reference frame of plane-parallel atmosphere. The z- 
axis is parallel to the unit vector n and the xy plane coincides 
with the plane of an atmosphere. The z'-axis (propagation di- 
rection of light) is the same as z-axis in Fig.[nand has an incli- 
nation 6 to the unit vector n and azimuth (p measured counter- 
clockwise from the x-axis: Q < 6 < n , Q < tp < In. The y angle 
has the same meaning as in Fig. [2 The magnetic field vector 
B is in the xz plane and has an inclination O to the unit vector: 
< Q < ;r. 



(40) 2.4. Optical depth scale 

In model atmosphere calculations with radiative transfer of po- 
larized radiation we can not use Rosseland depths tr as an op- 
tical depth scale because the extinction is not represented by a 
scalar but by a matrix. Moreover, this matrix depends on light 
propagation direction in the atmosphere, Eq. (jSJ. 

However, keeping in mind our assumption (0} that the con- 
tinuum is unpolarized it is convenient to introduce a continuum 
optical depth scale Tc, instead of using ts, by 



dTc — (kc + (Tc) dm. 



(45) 



We have used the monochromatic optical depth Tc calculated 
at 5000 A denoted here as rgooo, as a scale on which atmo- 
sphere model physical quantities are specified. The wavelength 
5000 A is usually used as a monochromatic optical depth in 
analysis of stellar atmospheres, and its choice is determined by 
the fact that the continuum opacity does not vary rapidly at this 
frequency and it is located around the maximum of the energy 
flow. Besides, the same depth scale is used in the LLModels 
code for solution of hydrostatic equation (^ulyak et al. 2 0041) . 

The model atmospheres calculations are carried out on a 
fixed T5000 grid using equal-space logrsooo scale subdivided 
into layers. 

2.5. Coordinate systems 

For a proper model atmosphere calculation in the presence of a 
magnetic field we need complete information about the orien- 
tation of the magnetic field vector B and the light propagation 
direction in the model reference frame. In Fig.|2]two cartesian 
references frames are represented. One of them (xyz) is related 
to the plane-parallel atmosphere, another, local reference frame 
(x'y'z') depends on the light propagation direction and the z'- 
axis coincides with this direction. The x' and y' axes are not 
shown in figure because their choice does not modify the main 
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equations used in this paper (Sect, \2.3l . Thus, three angles: 6 
and ip (inclination of the propagation direction to the unit vec- 
tor n, and its azimuth), and y (inclination of magnetic vector B 
to the propagation direction) in the model reference frame xyz 
fully represent the necessary orientation related information. 

The direction of the magnetic field vector in the model ref- 
erence frame xyz is set by the inclination angle D.. It is clear 
that an additional azimuthal angle is not necessary because of 
the integral nature of the flux and radiative equilibrium condi- 
tions, Eqs. (I42> and (I44> . That is why it is convenient to intro- 
duce the model reference frame xyz so that the x-axis is in the 
(B, n) plane. In this case the angle 7 is defined by formula 

cosy = cos cos O H- sin sin Q cosy?. (46) 

Consequently, one preset angle Q. and the spherical coordi- 
nate system with angles 6 and if completely describe the orien- 
tation of the magnetic field vector for computing models. 

3. Results of numerical application 

In order to test model atmospheres including polarized radia- 
tion transfer we calculated the same set of model atmospheres 
for A and B stars with the same modelling parameters as in 
Paper I. The main aim of this numerical application is to anal- 
yse effects of introduction of the polarized transfer equation in 
model atmosphere computation, because we suppose that the 
treatment of Zeeman splitting used in our previous work over- 
estimates the magnetic intensification. 

In Paper I we treated individual Zeeman components of 
the anomalous splitting pattern as independent lines by mod- 
ification of the original lines list and gf values for these lines. 
Splitting pattern incorporated into the line list supposes that the 
orientation of the magnetic field vector is characterized by the 
same angle y for every light propagation direction. That means 
that we neglected anisotropy, which is naturally arises in mag- 
netic atmospheres and is determined by the orientation of the 
magnetic field vector. In Paper I we assumed the angle y is al- 
ways equal to n/2 to minimize eff'ects of the polarized radiative 
transfer. That allowed us to use the transfer equation for non- 
polarized radiation and simplified the task. 

In other words, such an approach supposes that the mag- 
netic field orientation is assigned not with respect to the plane 
of atmosphere (or observer), but to the light propagation direc- 
tion. That is, the propagation matrix (|3 does not depend on 
angle 6, and angle y is fixed (Fig. |2j. On the one hand this is 
an unrealistic configuration of the magnetic field; on the other 
hand, for spherical stars, the angle dependent Eq. Q calculated 
in the plane-parallel approximation in an isotropic medium pro- 
duces variation of the intensity across visible stellar disk that is 
the same as the local angular variation. This applied magnetic 
field configuration was called "horizontal" in Paper I. 

Using the techniques detailed in Sect.|5]and taking into ac- 
count considerations outlined above to correspond to the previ- 
ous investigation, we have calculated model atmospheres with 
the eff^ective temperature = 8000K, llOOOK, 15 000K, 
surface gravity logg - 4.0, metalHcity [M/H] - 0.0, -1-0.5, 
-I- 1.0 (for non-solar metallicities the He abundance was de- 
creased to normalize the sum of all abundances) and magnetic 



field strength 0, 1, 5, 10, 20 and 40 kG, which covers a substan- 
tial part of the stellar parameter space occupied by the magnetic 
Ap and Bp stars. 

The model atmosphere calculations were carried out on the 
continuum optical depth grid spanning from +2 to -6.875 in 
IogT5ooo and subdivided into 72 layers. Convection was ne- 
glected as we discussed above. We also adopted zero microtur- 
bulent velocity in all calculations, because of direct inclusion of 
the magnetic intensification in the modelling of spectral lines 
and because the same approach was used in Paper I. 

At the start of the model computation for each of the 
^Rff, \M/H] pa irs studied, we have used a standard ATLAS9 
(lKuruczlll993ah model atmosphere for the spectral line prese- 
lection procedure and as an initial guess of the model structure. 
The line selection threshold was set to 1 % as previously. 

We used the iterative approach for the calculation of the 
magnetic model atmosphere grid presented here as described 
in Paper I. Each new model is calculated using a converged 
model with a weaker field strength (i.e. the OkG model is used 
as an initial guess for the 1 kG model, then the 5 kG model is 
calculated starting from the 1 kG model, etc.). 

The following wavelength ranges for the spectrum syn- 
thesis were used: from 500 A to 50 000 A for T^s = 8000 K, 
from 500 A to 30 000 A for T^g ^ II 000 K and from 100 A to 
30000 A for Teff - 15 000 K. For all models we used a wave- 
length step of 0. 1 A, which resulted in the total number of fre- 
quency poi nts in the ran^e betw een 295 000 and 495 000. It 
was shown dShuIvak et a l.' 2004) that a wavelength step 0.1 A 
is quite sufficient for model atmosphere calculations with the 
fundamental parameters used in this paper. 

In the following sections we describe preparation of the line 
lists, and we study the effects of the magnetic line blanketing 
on the common photometric and spectroscopic observables: the 
spectral energy distribution, photometric colors and profiles of 
the hydrogen Balmer lines. We do not consider model structure 
(i.e. temperature and pressure distribution) because we have 
used an optical depth scale different from that which was used 
in Paper I. We generally believe that all main trends concern- 
ing changes of model structure due to backwarming outlined 
in Paper I apply to this study too. Thus, those who are inter- 
ested in model structure changes due to Zeeman effect should 
be referred to Paper I (Sect. 3.1). 

3.1. Line lists and preselection procedure 

The technique for line list preparation is described in Paper I, 
and we used the same one. The magnetic line blanketing was 
taken into account for all spectral lines except the hydrogen 
lines according to the individual anomalous Zeeman split- 
ting patter n. The initial line list was extracted from VALD 
jKupka et al. 1999 ) including lines originating from predicted 
levels. The total number of spectral lines was more than 22.5 
million for the spectral range between 50 and 100 000 A. For 
0.55% of these lines information about Lande factors was ab- 
sent. Using the term designation represented in the VALD and 
assuming the LS coupling approximation, we computed Lande 
factors for the Unes of light elements, from He to Sc. That al- 
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Table 1. Description of the spectral line lists used in the line 
blanketing calculations for different effective temperature Teff 
and metallicity [MjH] = log(A^nietais/A^H) - log(A^inetais/A^H)0- 
Mines is the total number of spectral lines retained by the pres- 
election procedure and included in the line blanketing calcula- 
tion. 





[M/H] 


^lines 


8000 


0.0 


352 389 




0.5 


513 034 




1.0 


744 190 


11000 


0.0 


305 249 




0.5 


442 366 




1.0 


652 839 


15 000 


0.0 


361256 




0.5 


533 787 




1.0 


798 154 



lowed us to reduce the number of lines with unknown Lande 
factors (or without a proper term designation) to 0.26%, for 
which we assumed a classical Zeeman triplet splitting pattern 
with effective Lande factor geft - 1-2. The resulting line list 
was converted to special binary format file accepted by the 
LLMoDELS code. 

The resulting line list was used for the preselection pro- 
cedure in the LLModels code using the selection threshold 
1% that means that the code selects spectral lines for which 
{y/o^v > 1%, where (y and ay are the continuum and line ab- 
sorption coefficients at the frequency v. During the preselection 
procedure we have used a standard ATLAS9 ('Kurucz'lQgSa") 
model atmospheres for each of the Tefr, [M/H] pairs studied, 
assuming a non-magnetic case. 

Preselection allowed us to decrease the number of spec- 
tral lines involved in the line blanketing calculation to about 
300 000 - 800 000, depending on the model atmosphere param- 
eters. The number of preselected lines is presented in Tabled 

3.2. Energy distribution 

The emergent spectral energy distributions Fy (see Eq. i4-4l ) 
for the three values of effective temperature and for metallic- 
ity [M/H] = +0.5 and [M/H] = +1.0 are presented in Fig.|3] 
The theoretical energy distributions for magnetic stars are com- 
pared with the calculation for the reference non-magnetic mod- 
els. 

We found the same changes to the energy distributions 
compared to non-magnetic models as in Paper I. The first is 
a flux deficiency in the ultraviolet spectral region coupled with 
a flux excess in the visual. The magnitude of the UV deficiency 
is a moderate function of the effective temperature and clearly 
increases with the magnetic field strength and metallicity. 

The second related effect is that the presence of a magnetic 
field changes the stellar flux distribution in opposite directions 
in the visual and UV regions. At short wavelengths the mag- 
netic star appears to be cooler in comparison to a non-magnetic 
object with the same fundamental parameters, whereas in the 
visual magnetic star mimics a hotter normal star. The "null 
wavelength" where flux is unchanged progressively shifts to 



bluer wavelengths as the stellar effective temperature increases. 
The results of our modelling suggest that the "null wavelength" 
falls approximately at 3800 A for T^g = 8000 K, at 2900 A for 
Teft = 1 1 000 K and at 2600 A for T^k = 15 000 K, and depends 
slightly on metallicity and field strength. The same "null wave- 
lengths" were found in Paper I. 

We found that the model atmospheres with magnetic line 
blanketing produce fluxes that have a maximum deficiency in 
the UV region of about 0.1-0.2 mag for B - 10 kG and about 
0.2-0.3 mag for B = 40 kG, which is about half the effect found 
in Paper I. 

The third effect is that the theoretical flux distribution of a 
magnetic star exhibits a depression in the 5200 A region. This 
well-known depression is of interest because it is f requently 
obser ved in the spectra of peculiar stars and is used jMaitzenl 
to detect magnetic stars photometrically. In our theoret- 
ical calculations the 5200 A depression is prominent at lower 
Teff but becomes rather small for hotter models. The magni- 
tude of the depression increases with the magnetic field inten- 
sity and metal content of the stellar atmosphere. 

3.3. Colors 

We studied the influence of magnetic line blanketing on the 
photometric colors in the Stromgren uvbyWjS, Geneva and Aa 
systems. 

All colors were calculated using modified computer codes 
by iKurucz ( 1993a), which take into account transmission 
curves of individual photometric filters, mirror reflectivity and 
a photomultiplier response function. Our synthetic colors are 
computed from the energy distributions sampled every 0. 1 A, 
so integration errors are expected to be negligible. 

Having introduced a more accurate treatment of Zeeman 
splitting on line blanketing, it is particulary interesting to 
study specific photometric indices whose aim is to identify 
CP stars. Three such indices have been identified. The pecu- 
liar parameter Aa was introduced by Majtzen ( 1976 ) to mea- 
sure the strength of the 5200 A depression and is based on 
narrow-band filters. The A{V\ - G) parameter of the Geneva 
medium-band system can b e used as a peculiarity parameter 
for Ap stars llHau ck[|l274j). The Z parameter, introduced by 
ICramer & Maeder' ('l979'l is linear combination of the Geneva 
photometric indices and was suggested as an indicator of chem- 
ical peculiarity. All of these indices are sensitive to the depres- 
sion around 5200 A and to the strength of the surface magnetic 
field. 

The Aa system wa s specifically d esigned for 5200 A de- 
pression measurement (lMaitzenll976h . while the central wave- 
lengths of the Vi and G filters of the A(V\ - G) system are lo- 
cated approximately at 5300 A and 5800 A respectively, so the 
region of sensitivity of the A{Vi - G) system is shifted to the 
red in comparison to Aa system. The Z index is mainly sensitive 
at Vi filter ( Cra mer & Maeder 1980), so it is sHghtly shifted to 
the red, too. 

Here we compare our results with those derived in Paper I, 
relying on fact that the behaviour of the photometric indices is 
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Fig. 3. Synthetic energy distributions from UV to near IR region for effective temperatures Teff = 8000 K, 1 1 000 K, 15 000 K and 
metallicities [M/H] = +0.5, [M/H] = +1.0. To enable visual comparison of the different curves the original LLModels stellar 
fluxes were convolved with a Gaussian profile with FWHM = 15 A. The insets show energy distributions in the 5200 A region. 



closely related to the flux redistribution between the visual and 
UV regions and the presence of the flux depression at 5200 A. 

In order to save space for detailed illustrations of pecu- 
liarity indices we do not show general graphs for Stromgren 
uvbyH/3 and the Geneva photometric systems, because their are 
not particularly interesting. However, we note some common 
features for all photometric indices. Relations between their 
changes and the intensity of magnetic field depends strongly 
on the stellar effective temperature. For low Teff photometric 
changes are very pronounced, whereas for hotter magnetic stars 
modification of the photometric observables is fairly small (ex- 
cept ci). All photometric indices become more sensitive to the 
magnetic effects with increasing metallicity. These main out- 
comes coincide with results obtained in Paper I, except that the 



sensitivity of all non-peculiar photometric indices is about 30- 
40% less than those presented in the earlier paper. 

The ci of the Stromgren system measures the Balmer dis- 
continuity, which is why the shift of the "null wavelength", 
which leads to changes of the Balmer jump amplitude, affects 
the ci index. The relation between ci and B demonstrates the 
same, almost linear, behaviour as in Paper I but the overall ef- 
fect is approximately half as large, which reflects weaker flux 
redistribution (see Sect. 13. 2> . 

The results of our calculations of the synthetic peculiar a, 
Vi - G and Z photometric parameters are summarized in Fig.|4] 
and show their sensitivity to the metal abundance and magnetic 
field strength. Similarly to the effects observed for other pho- 
tometric indicators, these peculiar parameters are most influ- 
enced by magnetic field at lower T^s. In comparison to the 
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Fig. 4. The magnitude of the a, Vi - G and Z photometric indices as a function of the magnetic field strength for Tefr = 8000, 
11000 and 15 000 K. Different curves show calculations with metal abundances [M/H] = 0.0 (open circles), [M/H] = +0.5 
(triangles) and [M/H] = +1.0 (filled circles). Open squares illustrate the effect of increasing Cr overabundance to [Cr] - +2.0 in 
the Teff = 8000 K, [M/H] = +1.0 models. All curves are shifted so that for [M/H] - 0.0 and zero magnetic field in each frame 
the a,V\ - G and Z indices equal zero. 



results of Paper I the sensitivity of these indices to the mag- 
netic field strength appears to be between 10-50% less as we 
decrease T^s from 15 000 K to 8000 K. 

The Afl system is most sensitive to the magnetic field 
strength, except for model atmospheres with T^s - 8000 K 
where the I^{V\ - G) demonstrates the best sensitivity. 
Furthermore, it seems that the Aa system is not very sensitive 
to metallicity (although it does increase with increasing B) for 
low Teff. While for non-magnetic case the small difference be- 
twee n a i ndices for metallicities other then solar is normal fea- 
ture jKupka et al. 2003. Fig. 4), the small variation of a with 
[M/H] for strong magnetic field is not expected (in Paper I we 
found strong dependence of the a-index on metallicity for all 
T^eff). We suppose that this phenomena appears due to some sat- 
uration effect of enhanced line blanketing in the narrow region 
around 5200 A, which is mainly dominated by the Fe I and low 
excitation Fe II lines for low effective temperature. In contrast 
to the narrow-band Aa-system the sensitivity of medium-band 
parameters V\ - G and Z does not suffer much from the spe- 
cific properties of spectral feature in the naiTow interval around 
A 5200 A. This state is well illustrated in flux units (Fig.|5jl. One 
can see that flux depression in the region longward of 5200 A 
is clearly affected by additional abundance and magnetic in- 



tensification, while the narrow region centered at 5200 A is not 
altered much. 

To test the influence of an anomalous concentration of an 
individual species on the photometric peculiarity parameters, 
and especially the Aa value for low effective temperature, we 
increased the Cr overabundance by a factor 100 relative to 
the sun and calculated additional model atmospheres. As il- 
lustrated in Fig. 13 and Fig.lsjthis results in 0.012-0.015 mag 
growth oftheAa index for Teff = 8000 K, [M/H] = +1.0 model 
and field strength 5-20 kG. For A(yi - G) the increase is 
smaller and amounts to 8-13 mmag while, the Z-index is al- 
most unaffected (1-3 mmag). This behaviour indicates that 
observable anomalous Aa values for low effective temper- 
ature may arise not only from the influence of the mag- 
netic field but also from specific anomalous abundances as 
well. Nevertheless, even with overabundant Cr the Aa system 
demonstrates the least sensitivity to the magnetic field strength, 
with only 8 mmag change from 5 to 20 kG, while the value for 
Z is 13 mmag and 20 mmag for A{V\ - G). 

Finally, none of the photometric peculiarity indicators (ex- 
cept a for Teff = 8000 K) show a linear trend over the whole 
range of the considered magnetic field strength. Saturation ef- 
fects for their changes with the magnetic field intensity ap- 
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Fig. 5. Energy distribution around 5200 A for effective temper- 
ature Teff - 8000 K, B - 20 kG and two values of metallicity 
[M/H] - +0.5 and [M/H] = +1.0; for second metallicity the 
behaviour for overabundant [Cr] - +2.0 is also shown. The top 
frame represents the comparison of fluxes from present study, 
the bottom frame shows the respective data from Paper I. The 
fluxes are convolved with a Gaussian profile assuming reso- 
lution R = 150. The curves of photometric filters (dash-dotted 
lines) are presented for comparison. The gi, g2 and y are pro- 
files of filters of the Aa system, and Vi and G of the Geneva 
system. 



pear for B > lOkG. The Aa index for low eff'ective tempera- 
ture (Teff = 8000 K) demonstrates almost linear dependency on 
B from 5 to 40 kG. 



3.4. Hydrogen line profiles 

Profiles of the Ha, H/3 and Hy hydrogen Balmer lines 
were calculated using the Synth spectrum synthesis program 
dPiskunovlllQQl . The comparison between the magnetic and 
non-magnetic profiles of the li/3 line calculated for models with 
[M/H] - +1.0 metallicity is presented in Fig.|6l 

We found that the introduction of the polarized radiative 
transfer to the magnetic line blanketing calculation do not have 
a very strong influence on the Balmer line profiles. For proper 
further discussions we note that hydrogen line profiles calcu- 
lated in this work for B - OkG are almost the same as those 
calculated in Paper 1 (within < 1% error bar). 

For all Balmer lines considered, the maximum change 
relative to the non-magnetic model does not exceed 1% for 
B < 10 kG, which is true of the surface field strength for the 
great majority of magnetic CP stars, while the same quantity 
presented in Paper 1 (Sect. 3.4) is less then 2%. 

For metallicity [M/H] = 0.0 the maximum change among 
aU Balmer profiles relative to the non-magnetic models reaches 



0.5% of the continuum level for B < lOkG and 1.0% for 
B = 40 kG. 

For the overabundant composition ([M/H] - +1.0) the dif- 
ference between hydrogen profiles for magnetic and non- 
magnetic models may be as large as 1.0% for B < lOkG and 
1.8-2.5% for B = 40 kG. 

For the li/3 line the maximum change amounts to about 2% 
of the continuum level for B = 40 kG. The maximum devia- 
tion of 2.5% found in the whole grid of models occurs for Ha 
(Feff = 11 000 K, [M/H] = +1.0, and B = 40 kG), which ex- 
hibited the most sensitivity to magnetic field effects in Paper I, 
too. 

3.5. Stellar atmospheric parameters and bolometric 
correction 

The modification of the energy distribution due to magnetic 
line blanketing may change photometric determination of the 
stellar atmospheric parameters. To verify its possible influ- 
ence o n Teff and log g determination we applied the TempLogg 
code l| R9g?r^ l^ ??-'^^ w hich implements the calibration by 
"Napi wotski et aLldiggsl) to synthetic uvbyH/3 colors (Sect.|33}. 

As in Paper I we examined the error introduced by the mag- 
netic field on the determination of the stellar atmosphere pa- 
rameters. For each set of the eff'ective temperature and metallic- 
ity we calculate differences between the values of Tf^g and log g 
for non-magnetic and for the Z? = 40 kG models, then the max- 
imum deviations of T^ff and log g among different metallicities 
are chosen. We found that for Tss = 8000 K the maximum error 
for Teff and log g introduced by 40 kG magnetic field is-128K 
and -0.14dex respectively. For the Tgff = 11 000 K models the 
temperature discrepancy is +8 1 K and for log g is -0.29 dex. 
For the hottest models (Teff - 15 000 K) the Tgg difference is 
+502 K, and for \ogg is -0.12 dex. 

At the same time the overabundant models demonstrate in 
most cases higher values of the effective temperature and \ogg, 
but the differences in comparison to the solar composition mod- 
els do not exceed 104 K and 0.24 dex for T^ft = 8000 K, 106 K 
and 0.13dex for T^s = 11 000 K, and 538K and 0.15dex for 
Teff = 15 000 K. 

Finally, the differences between maximal and minimal 
values of the T^ff and \ogg for all values of the magnetic 
field strength and metallicity for T^ff - 8000 K are 152 K 
and 0.32 dex, respectively, for Tgff = 11 000 K are 151 K and 
0.29 dex, and for Teff = 15 000 K are 868 K and 0.16 dex. 
The same values but for B < lOkG are smaller and equal to 
104 K and 0.25 dex for T^g = 8000 K, 87 K and 0.09 dex for 
Teff = 11000 K, and 5 13 K and 0.15 dex for Teff = 15000 K. 

We confirmed as in Paper I that even in the case of extreme 
magnetic field strengths, the bias in determination of the funda- 
mental parameters using Stromgren photometry does not result 
in a systematic error in Teff and \ogg beyond the usual error 
bars assigned to the photometrically determined stellar param- 
eters. 

We have not conducted the same estimate of possible in- 
fluence on Tetf and log g using Geneva photometric colors be- 
cause in Paper I was shown that this system is less robust for 
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determination of the atmospheric parameters of magnetic stars 
compared to the calibration in terms of the MvfcyHyS colors. 

Due to redistribution of the stellar flux from UV to visual, 
magnetic CP stars appear brighter in V compared to normal 
stars with the same fundamental parameters. Using the theo- 
retical flux distributions we studied the effect of magnetic field 
and metallicity on the bolometric correction (BC). For all stel- 
lar models this parameter changes relative to the non-magnetic 
case by up to 0.029-0.053 mag for B - 40 kG and by no more 
than 0.018 mag for B < lOkG. However, in the absence of the 
field, an increase of metallicity changes BC. For instance, for 
metallicity [M/H] = +1.0 the change is 0.084-0. 125 mag. 

Generally, there no substantial differences in determination 
of atmospheric parameters and bolometric correction in com- 
parison to the results obtained in Paper I. 

4. Conclusions 

This paper continues our study of magnetic line blanketing in 
stellar atmospheres and its relation to the observable quanti- 
ties of magnetic CP stars. We presented a technique of model 
atmosphere calculation, introducing the transfer equation for 
polarized radiation. We have used a modified version of the 
model atmosphere code LLModels, which implements a direct 
treatment of the line opacity. 

Using this new technique, we computed a grid of model at- 
mospheres of peculiar A and B stars which are similar to those 
computed in Paper 1 in order to analyse effects of the introduc- 
tion of the polarized transfer equation. We have used this grid 
to analyze the behaviour of observed characteristics for differ- 
ent magnetic fields, abundances and effective temperature. To 
insure accurate modelling, we have used up-to-date atomic line 
data from the VALD database. This approach represents the 
most accurate treatment of the Zeeman effect in model atmo- 
sphere calculation used until now. 

We found that the model atmospheres with magnetic line 
blanketing produce fluxes that are deficient in the UV region, 
and the presence of a magnetic field leads to flux redistribution 
from the UV to visual region. This property of the theoretical 
models is in agreement with observations of CP stars and with 
results presented in Paper I, but the overall effect is approxi- 
mately twice smaller compared to the results of Paper I. 

We also investigate the most prominent feature of the 
CP stars in the visual region, the flux depression centered at 
5200 A. Our numerical results show that the magnitude of 
the feature near 5200 A depends strongly on magnetic field 
strength, T^g and metallicity. The depression grows with in- 
creasing magnetic field strength and metalUcity. However, for 
higher effective temperatures the magnitude of the depression 
for overabundant models becomes smaller, while the sensitivity 
to the magnetic field strength appears almost negligible. This 
outcome coincides with the conclusion found in the previous 
study. 

In order to study the flux depression around 5200 A more 
carefully we examined the influence of the magnetic field on 
the photometry observables, in particular, photometric pecu- 
liarity parameters Aa, Z and A{Vi - G), that are sensitive in 
this spectral region. We found the same trends as in Paper I, 



i.e. changes of peculiar parameters due to the influence of the 
magnetic field are noticeable for low effective temperatures, 
whereas for hotter stars sensitivity to the magnetic field is re- 
duced considerably, and the relation between B and the pecu- 
liarity index depends on metallicity. 

The magnetic modification of these indices is less by 30- 
50% (from higher to lower effective temperature) compared to 
those presented in Paper I. The Aa system is still most sensi- 
tive to the magnetic field (the maximum change is 22mmag 
for the range of B from to 40 kG and [M/H] = +1.0) for 
higher = 11 OOOK or 15000K while for T^b = 8000 K 
the A(Vi - G) system is preferable (the maximum change is 
38 mmag vs. 20 mmag for Aa for field strength running from 
to 40 kG and [M/ff] = +1 .0). The Aa system seems is not able 
clearly distinguish CP stars with metallicity other then solar for 
low effective temperature. Moreover, for low effective temper- 
afiire (TefF = 8000 K) and B = lOkG (which is typical for flie 
most strongly magnetic CP stars) it exhibits overall changes 
of only 9-13 mmag, while A{Vi - G) and Z for show overall 
changes of about 8^5 mmag. 

Our numerical tests show that Cr appears to be important 

contributor to the narrow wavelength interval around 5200 A 
where the Aa system is quite sensitive, while the A{V\ - G) 
and especiaUy Z system are almost unaffected by Cr overabun- 
dance. But even with overabundant Cr the sensitivity of the Aa 
system on magnetic field strength is the worst among other pe- 
culiar indices for low effective temperature. 

The relation between values of peculiar parameters and the 
magnetic field strength is not Unear due to the saturation ef- 
fects for stronger fields {B > 10 kG), except Aa which appears 
to have approximately linear dependance on magnetic field 
strength from 5 to 40 kG for Teff = 8000 K. 

Finally, we investigated the question of stellar atmosphere 
fundamentals parameters determination based on calculated 
synthetic photometric colors of the Stromgren system. We 
found that the model atmosphere parameters derived using the 
photometric calibrations for normal stars are not far from their 
true values and are within the usual error bars. We also showed 
that the maximum change of profiles of the hydrogen B aimer 
lines in comparison to the non-magnetic case does not exceed 
2.5% of the continuum level for the magnetic field strength 
40 kG. 

In our further research we plan to build model atmospheres 
of CP stars introducing a grid of 1 -D models on the stellar sur- 
face with individual values and orientations of magnetic field 
vector using technique described in this paper. This wiU allow 
us to study anisotropy effects induced by the magnetic field on 
line blanketing. 
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strength and Teff. The lower panels show normalized profiles, whereas the upper plots illustrate the difference between Hfi 
calculated for the models with substantial magnetic line blanketing and the reference non-magnetic model. 
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